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Abstract 

We present an investigation of the sinh-cosh (SC) interaction model with 

twisted boundary conditions. We argue that, when unlike particles repel, 

the SC model may be usefully viewed as a Heisenberg-Ising fluid with moving 

Heisenberg-Ising spins. We derive the Luttinger liquid relation for the stiffness 

and the susceptibility, both from conformal arguments, and directly from the 

integral equations. Finally, we investigate the opening and closing of the 

ground state gaps for both SC and Heisenberg-Ising models, as the interaction 

strength is varied. 
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In Ref. and 0, to be called I and II in the following, we solved the integrable one- 
dimensional (ID) SC-model defined by the Hamiltonian 

1 d 2 

H =~0 E ^2+ E Vjkfe-Xk). (1) 

1<3<N j l<j<k<N 

The pair potential is given as 
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s > -1, (2) 



and the quantum number a = ±1 distinguishes the two kinds of particles. We may usefully 
think of it as either representing charge or spin. For values of the interaction strength s in 
the range — 1 < s < 0, the system exhibits two gap less excitation branches with different 
Fermi velocities as does the repulsive ID Hubbard model and thus may be classified as a 
typical two-component ID Luttinger liquid [|J. The asymptotic behavior of the correlation 
functions is given by finite-size arguments of conformal field theory. A Wiener-Hopf type 
calculation shows that the spin-spin part of the dressed charge matrix is essentially 
identical to the dressed charge scalar in the Heisenberg-Ising (H-I) model ||. 

In this Letter, we will further explore the connection of the SC model with the H-I model 
by examining the response of the system to a flux <3>. The addition of a flux is compatible 
with integrability and allows the study of the transport properties by an adiabatic variation 
of $. For the H-I model, this has already been done ||[7| for the interaction strength 
range — 1 < A < 1. We will show that the spin degrees of freedom of the SC model for 
> s > — 1 may be usefully viewed as a H-I model with moving H-I spins. The presence of 
the translational degrees of freedom will simply renormalize the spin-spin coupling. 

We thus restrict ourselves in what follows to the unbound case — 1 < s < 0, such 
that there are two gapless excitations corresponding to a particle-hole and a two spin- 
wave continuum with excitation velocities v and v s , respectively. Let us then modify the 
Bethe ansatz equations of Eq.(II.7) by threading them with a flux $. We have two coupled 
equations for N particles with pseudo-momenta k = (ki, . . . , k^) and M spin waves with 
rapidities A = (Ai, . . . , Am) on a ring of length L. The energy of a given state is E(k) = 



2 Sj=i ^| an d the total momentum is -P(k) = J2f=xkj- Boosting the system by $ will 

accelerate the two kinds of particles in opposite directions due to the two components being 

of equal but opposite charge. Therefore, we have no center-of-mass motion and P = 0. The 

energy of a given state will change as a function of $, and the energy shift of the ground 

state may be written as A£ ($) = E ($)-E (0) = D$ 2 /2L + 0{^), where D is called the 

stiffness constant and can be specified by perturbation arguments for $ up to 7r ||. Note 

that since we do not have any center-of-mass motion, we can call D either spin or charge 

stiffness depending on what interpretation of a we adopt. We choose the spin language 

for comparison with the H-I model. However, the charge interpretation is probably more 

natural to describe transport properties. We furthermore caution the reader that the term 

charge stiffness has been previously used in lattice models to describe center-of-mass motion. 

The twisted Bethe ansatz equations are given by 

M M N 
- Lk 3 = 2nlj(kj) - — $ + ]T 0o,-i(*i - Aa) + E M*i - k d, (3a) 

iV 0=1 1=1 

M N 

= 27rJ a (A a ) + $ + - A 6 ) + XX-i(A a - kj). (3b) 

6=1 J'=l 

The two-body phase shifts for particle-particle, particle-spin wave and spin wave-spin wave 
scattering, 9 0t0 (k), 9 _i(k) and 6Li 5 _i(fc) respectively, have been given in I. The particle 
quantum numbers Ij and the spin-wave quantum numbers J a are integers or half-odd integers 
depending on the parities of N, M as well as on the particle statistics For simplicity, we 
use bosonic selection rules, although a purely fermionic or a mixed bose-fermi system may 
be studied along similar lines. In the ground state of the bosonic system, we have 

(N-l) (iV-3) (N-l) 
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(M-l) (M-3) (M-l) 
J 1 ,J 2 ,...,J M = , , , [4) 

for both N and M even. 

We start with some general considerations. Let us denote by E{ijy($) the energy of a 
state specified by the $ = set of quantum numbers {I, J}. We then adiabatically turn 



on the flux until we return to our initial state. The energy will also return to its initial 
value, although, it may return sooner; so, the period of the wave function will be an integer 
multiple of the period of the energy. We can define a topological winding number n to be 
the number of times the flux $ increases by 2ir before the state returnes to its initial value. 
As Sutherland and Shastry have shown, the ground state winding number of the H-I model 
with S z = in the parameter range — 1 < — cos(//) = A < 1 is 2, implying charge carriers 
with half the quantum of charge, except at isolated points A = cos(tc/Q), where M > Q > 2 
is an integer. In particular, at A = 0, the free particle wave function has periodicity 2ttNhi, 
where Nhi is the number of H-I sites, implying free acceleration in the thermodynamic limit. 

We now note the following important fact: Choosing p, = — its, the spin wave-spin wave 
phase shift 0-1,-1 is identical to the spin-spin phase shift in the H-I model, and we may 
rewrite the equation for the rapidities as 

N M 

Nd ,-i(\a, = JV £ 0o,-i(A a - k v ti/N = 2nJ a (X a ) + $ + £ 6>_i _i(A a - A 6 , //). (5) 

3=1 6=1 

which nearly is identical to the Bethe Ansatz equation of the H-I model, as can be readily 
seen when we use the standard transformation for the H-I momenta p = f(a,p). We then 
merely have to identify a = 7rA. The sole effect of the pseudo-momenta k is an averaging 
on the left hand side. Let us now restrict ourselves in what follows to the neutral (spin 
zero) sector such that we have M particles with a = — 1 and M particles with a = +1 
for a total of iV = 2M. Then, a discussion of the behavior of the rapidities A for varying 
$ exactly mimics the discussion of the H-I momenta p in Ref. at S z = 0: As long as 
|<3>| < 2tt(s + 1), all A's stay on the real axis. At $ = 2ir(s + 1), the largest root Am goes to 
infinity. For $ increasing beyond this point, \m will reappear from infinity as iit + 71 until 
exactly at $ = 2ir, Am = iir (7! = 0) and the remaining M — 1 rapidities have redistributed 
themselves symmetrically around on the real axis. However, as mentioned above, this 
behavior is different at the threshold values s — (1 — Q)/Q. The momenta k are always real 
and distributed about the origin. Eq.(|3I]) simplifies at $ = 2ti(s + 1) (and thus Am = 00) 
and is in fact just the equation for M — 1 rapidities in the ground state. So as in Ref. 
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using simple thermodynamical arguments, we may write 

AE (2n(s + 1)) = E (N, M - 1) - E (N, M) = 1/2L ■ X ~\ (6) 

where % * s the susceptibility. Comparing with the definition of the stiffness constant D, we 
find D = x~747r 2 (s + l) 2 . 

On the other hand, we can read off the finite-size energy corrections for the SC model, 
and then finite-size arguments of conformal field theory give an expression for AEq(2tt(s + 1)) 
in terms of the conformal weights, the dressed charge matrix H and the spin wave velocity 
v s . The neutral sector dressed charge matrix is given in Eq. (11.35) and thus we have = 
2irv s (s + 1). We may therefore express the stiffness D in terms of the spin wave velocity as 

D = v s /2n(s + l). (7) 

We emphasize that this formula for D is true also for a system of purely fermionic particles. 
Shastry and Sutherland || have given an exact formula for the stiffness constant in the H-I 
model, by using the known expression for the H-I model spin wave velocity v s = 7rsin(//)//i 
||. No such expression is known for the SC model and we can only give v s as 

1 J B B e-^ 2s e'(k)dk 
Vs 2n f* B e ~^/ 2s p(k)dk ' [ ' 

Here, we use the definitions of II, Section II. However, written in terms of spin velocities the 
stiffness formulas are identical and only the values of the respective spin wave velocities are 
different. Thus the presence of the translational degrees of freedom in the SC model simply 
renormalizes the spin-wave velocity. 

We have iterated the Bethe Ansatz equations (^) in the neutral sector for reasonably 
large systems and density d = N/L = 1/2 as a function of <3>. By our correspondence 
between the H-I model, and the spin wave part of the SC model, we expect free spin waves 
at s — — 1/2. In the thermodynamic limit, we would thus expect the periodicity of the 
ground state energy to be infinite. For a finite system, this will be reduced to a periodicity 
that scales with the system size. For the SC model we have indeed found that at s = —1/2 
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the periodicity of the ground state energy is given as 2irN. We may then speak of s — > — 1 + 
as the ferromagnetic critical point and s — > 0~ as the antiferromagnetic critical point of the 
SC model. In Fig.(|I]) we show the full spectrum of low-lying states with zero momentum 
at s = — 1/2 for L = 12, iV = 6 and M = 3. The ground state curve is emphasized and its 
periodicity is 6 • 2n. 

Note that at $ = 2ir there is a level crossing between ground state and first exited state 
in Fig. ([!]). When the interaction strength changes from s = —1/2, immediately a gap opens 
between the ground state and first exited state. Just as in the H-I model the periodicity is 
reduced to Air. Note that a perturbation theory argument can not describe this behavior. 
Fig.(Q) shows the behavior of the ground state energy variation L[l — E($)/E(2ir)] for 
s = —1/3 near $ = 2n for different lattice sizes. The rounding is well pronounced and does 
not vanish as we increase the size. 

Thus the behavior of the low-lying states in the SC and H-I models is qualitatively the 
same, up to the renormalization of quantities such as the spin wave velocity v s . Let us briefly 
describe the behavior of the gaps in the H-I model, keeping in mind the correspondence 
\i = —us. Increasing /i beyond 7r/2 (A = 0), we see that the gap continues to widen up to 
a maximum value at \x ~ 77r/12 (A ~ 0.26). It then closes up again exactly at /i = 27r/3 
(A = 1/2). As has been noted before, this value of \i coincides with the appearance of 
a Q = 3 string ||. Further increase of /J again opens, and then closes the gap at the 
threshold for the next-longer Q = 4 string. This behavior continues, and the threshold 
values accumulate as /x — > n (A — > 1). In Fig.(|3|), we show the ground state and the first 
exited state of the H-I model on a ring of Nhi = 12. Note that due to the finite size of 
the ring, we can only observe strings up to length Q = 6. We will present a more detailed 
finite-size study of the behavior of the gaps in H-I and SC model in another publication. 
We only mention that for fixed fi the gap scales with the system size as a negative power of 
Nhi, with variable exponent depending on the coupling constant \i. 

The stiffness constant D is the curvature of the ground state energy Eq(Q) as a function 
of $. In Fig.(U), we show D for systems of 12, 24 and 32 lattice sites. We also show the 
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behavior of D as given by Eq.(|7]). As s — > _ , the spin wave velocity approaches the velocity 
of a non- interacting single- component model, i.e. v s — > nd/2 [Q]. Thus D approaches the 
non-zero value 1/8 which is compatible with the result of Ref. ||. Furthermore, the SC 
model exhibits a gap for s > and so D is zero. Thus D exhibits a jump discontinuity at 
s = just as in the H-I model for A = — 1. 

Note that Eq.flT]) may also be written as D\ 1 = v 2 s . This is nothing but the Luttinger 
relation for the spin wave excitations ||. Let us briefly explain how to derive this formula 
without using arguments of conformal field theory. In the thermodynamic limit, we convert 
Eq.(^|) into a set of coupled integral equations as in Eq.(II.lO). Here p{k) and <r(A) are 
the distribution functions of particles and down-spins, respectively. The density d and the 
magnetization M. are then given parametrically in terms of the integral limits B and C. We 
now use an iteration scheme, i.e., first, for B = oo and M. = 0, i.e. at half-filling and zero 
magnetization, we calculate p{k) = po(k). We then use this po(k) in the equation for cr(A) 
with C finite and B nearly oo. Finally, we use this cr(A) to calculate p(k) and thus the effect 
on the momenta and the energy. Since we are only interested in the leading order correction 
terms, we may stop. Then the corrections to the energy are 



AE _ 1 

~zT = 2 



v. 



(|) 2 + [27r(l + S )] 2 ^ 2 



(9) 



4tt(1 + s 

where M. = \d(l — 2M/N). A complete account of this calculation can be found in Ref. 0. 

The derivation of the Luttinger relation uses integral equations and as such is valid 
in the thermodynamical limit. Most of the other results given above have been derived 
using Eq.(^[). These equations, however, have been derived by the asymptotic Bethe Ansatz 



(AsymBA). This method is only correct in the thermodynamic limit ||10|| . Thus all our finite- 
size results should exhibit correction terms. From the hyperbolic form of the pair potential 
(0), we may expect these corrections to be exponentially small in L. Indeed, a log- log plot 
of the ground state energy versus L at fixed interaction strength shows a simple power law 
behavior already for L > 6. Thus the L — > oo behavior of the finite-size Bethe- Ansatz 



equations for the SC model does not seem to differ in any significant respect from usual 
finite size behavior for short ranged models. This further supports our use of the AsymBA 
in the present study. 

In conclusion, we have shown that the SC model exhibits all the rich structure of the H-I 
model for — f < s < 0. In particular, there is a Luttinger relation for the spin waves just as 
in the H-I model, that can be derived from (i) conformal arguments, (ii) an exact calculation 
in the thermodynamic limit and (iii) is furthermore supported by numerical results for finite 
systems. Thus this yields credibility to both the conformal and the Luttinger approach in 
models solved by the AsymBA. Finally, we have reported an interesting behavior of the gaps 
in H-I and SC models. 

R.A.R. gratefully acknowledges partial support by the Germanistic Society of America. 
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FIGURES 

FIG. 1. The low-lying states for the SC model at L = 12, N = 6 and M = 3. The bold curve 
corresponds to the ground state and the winding number is n = 6 = N. Note the various level 
crossing in this free spin wave case, especially the crossing of ground state and first exited state at 
$ = 2vr. 

FIG. 2. Plot of the ground state energy variation L[l — E(&)/E(2n)] for the SC-model at 
s = -1/3 for L = 12, 20 and 28. 

FIG. 3. The charge stiffness D(s) for the SC model. The dashed curves correspond to L = 12, 
24 and 32 and converge to D(0) = 1/8 at s — > 0~. The solid curve comes from Eq.(^), which can 
be derived by conformal methods or from the Luttinger relation. (Note that as s — > 0~, the solid 
curve does not converge to 1/8. This is due to a buildup of numerical errors in the integration 
routine.) 

FIG. 4. Energy of the ground state and first exited state and their difference in the H-I model 
for Nhi = 12. Note the closing of the gap at A = cos(7r/Q) for Q = 2, 3, 4, 5. 
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